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made  to  determine  the  oxlsymmetrlc  stress  In  the  solid  Induced  by  a  sudden 
pressure  distribution  during  the  collapsing  of  a  bubble. 

An  existing  thln-fllm  analysis  for  determining  the  flow  faictors  used  In 
the  hydrodynamic  analysis  of  rough  surfaces  was  extended  to  Include  the  steady 
state  cavitation  effects.  Results  showed  that  the  shear  flow  factors  0*  cal¬ 
culated  using  a  model  considering  cavitation  are  dependent  on  ambient  pressure 
surface  pattern  parasieter  and  the  film  thickness  to  roughness  ratio,  h/a. 

The  model  predicts  a  lower  shear  flow  factor  when  cavitation  Is  con¬ 
sidered.  They  approach  asymptotically  to  the  non-cavltated  values  as  the 
aniblent  pressure  rises.  The  effect  of  cavitation  are  found  to  be  less  signi¬ 
ficant  for  h/a  >3or  for  h/o  >  0.75  and  more  significant  for  the  Intermediate 
values.  The  surfaces  with  transverse  roughness  resulted  In  more  cavltated 
regions  than  Isotropic  and  longitudinal  roughness. 
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ABSTRACT 


An  experimental  rig  consisting  of  a  hydrostatic  circular  step  pad  was 

constructed  to  simulate  the  cavitation  phenomenon  under  a  water  pressure  up  to 

2000  psl  comparable  to  those  In  the  submarine  stem>tube  seals.  Cavitation 

bubbles  were  observed  for  water  glycerine,  and  their  mixtures. 

Several  analytical  models  were  studied,  and  the  physical  parameters 

Influencing  bubble  growth  and  collapse  were  Identified.  An  analysis  was  also 

made  to  determine  the  axlsymmetrlc  stress  In  the  solid  Induced  by  a  sudden 

pressure  distribution  during  the  collapsing  of  a  bubble. 

An  existing  thln-fllm  analysis  for  determining  the  flow  factors  used  In 

the  hydrodyamlc  analysis  of  rough  surfaces  was  extended  to  Include  the  ste^'dy- 

* 

state  cavitation  effects.  Results  showed  that  the  shear  flow  factors  6  cal- 

8 

culated  using  a  model  considering  cavitation  are  dependent  on  ambient  presure, 
surface  pattern  parameter  and  the  film  thickness  to  roughness  ratio,  h/o  . 

The  model  predicts  a  lower  shear  flow  factor  when  cavitation  Is  con¬ 
sidered.  They  approach  asymptotically  to  the  non-cavltated  values  as  the 
ambient  pressure  rises.  The  effect  of  cavitation  are  found  to  be  less  signi¬ 
ficant  for  h/o  >  3  or  for  h/o  >  0.75  and  more  significant  for  the  Intermediate 
values.  The  surfaces  with  transverse  roughness  resulted  In  more  cavltated 
regions  than  Isotropic  and  longitudinal  roughness. 


I.  nmtoDucTioN 

In  high  pressure  face  seals,  such  as  the  submarine  stem-tube  seals,  It  Is 
often  found  from  the  prematurely  failed  seals  that  sections  of  the  seal  surface 
show  a  matted  appearance  similar  to  the  surfaces  damaged  by  cavitation  (1).  This 
mode  of  failure  appears  to  be  even  more  prevalent  as  both  the  pressure  differential 
and  the  seal  size  are  enlarged.  It  has  been  suggested  that  such  surface  damage 
may  be  caused  by  the  collapsing  of  cavltated  bubbles  associated  with  the  sudden 
drop  In  pressure  In  a  rapidly  diverging  channel.  In  order  to  understand  the 
basic  mechanism  of  cavitation  damage  In  a  narrow  gap  between  the  sealing  surfaces, 
a  research  program  was  Initiated  at  Northwestern  to  meet  the  following  objectives: 

a)  Design  and  construct  an  experiment  to  simulate  the  cavitation  phenomenon 
under  conditions  comparable  to  those  In  the  recent  high  pressure-stern- 
tube  seals. 

b)  To  formulate  an  analytical  model  to  explain  the  growth  and  collapse  of 
the  cavltatlng  bubbles  In  a  narrow  film. 

c)  To  employ  an  existing  thin- film  analysis  for  rough  surface  to  predict 
the  effect  of  cavitation  on  the  steady-state  performance  between 
sliding  sealing  surfaces. 

This  report  presents  the  results  obtained  during  the  first  phase  of  the 
above  described  tasks. 


II.  CAVITATION  IN  FACE  SEALS 

Mechanical  seals  consist  of  two  surfaces  In  relative  motion  and  In  intimate 
contact.  The  gap  between  the  parts  is  of  the  order  of  the  surface  irregularities 
(typically  a  fraction  of  a  micrometer).  The  formation  and  collapse  of  vapor 
bubble  in  the  fluid  that  fills  the  thin  sealing  gap  is  known  as  cavitation  in 
face  seals. 

Cavitation  in  external  flows;  e.g.,  around  propellers  or  moving  objects 
immersed  in  water.  Involves  vapor  bubbles  growing  in  the  liquid  at  low  pressure 
regions,  and  subsequently  collapsing  when  exposed  to  a  higher  pressure.  Experi¬ 
ments  have  shown  these  bubbles  start  as  nuclei  in  the  water  of  about  1  micrometer 
in  size  (2) .  Nuclei  of  this  size  in  a  seal  gap  must  grow  into  disk  shaped  vapor 
bubbles  because  of  the  constraining  surfaces.  The  fundamental  geometry  is 
cylindrical  rather  then  spherical  (Fig.  1  illustrates  this  difference) . 

Cavitation  damage  to  solid  surfaces  is  associated  with  vapor  bubbles 
collapsing  in  the  neighborhood  of  the  surface  (3).  In  order  to  explain  cavitation 
damage  in  seals,  it  is  necessary  to  establish  how  vapor  bubbles  grow  and  collapse 
in  the  gap  region.  A  bubble  must  grow  to  a  certain  size  before  its  collapse  can 
generate  enough  disturbance  to  cause  damage. 

Several  paths  can  produce  this  history  for  a  bubble.  The  most  direct  occurs 
when  the  bubble  remains  stationary  in  the  fluid  and  the  local  pressure  first 
decreases  and  then  increases.  Such  a  pressure  history  might  be  generated  by 
squeeze  film  action  if  the  gap  between  the  two  surfaces  periodically  changes 
(due  to  vibration  or  other  dynamical  effects) .  For  unconflned  flows  cavitation 
is  produced  in  this  manner  by  ultrasonic  generators.  Unsteady  pressure  fields 
such  as  these  are  also  produced  by  roughness  on  the  two  surfaces.  As  asperities 
pass  one  another,  the  neighboring  fluid  is  subjected  to  a  time  varying  pressure 
field. 
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Another  path  Is  through  convection,  wherein  a  bubble  Is  formed  in  a  low 
pressure  region  and  is  then  transported  by  the  fluid  to  a  region  idiere  the  pres¬ 
sure  is  higher.  A  typical  example  of  this  process  for  unconfined  flows  occurs 
on  submerged  bodies  where  vapor  bubbles  originate  at  low  pressure  shoulder  regions 
and  collapse  when  transported  to  a  reattachment  region  (4). 

An  important  question  in  this  case  is  how  a  steady  flow  produces  an  unsteady 
flow  of  bubbles.  Turbulence  plays  a  role  as  well  as  the  stability  properties 
of  the  interface  between  the  vapor  and  liquid.  The  incompatabillty  of  the  con¬ 
stant  pressure  vapor  filled  region  and  the  adjacent  liquid  flow  which  requires 
pressure  changes  along  the  boundary  is  also  important  (5).  This  latter  effect 
mi^t  also  be  Important  in  the  seal  problem,  whereas  turbulence  and  interfacial 
instability  are  unlikely  to  be  of  importance  because  of  the  characteristically 
low  Reynolds  number  for  seal  flows. 


III.  STEADY  VAPOR  BUBBLES  IN  SEAL  FLOW 


The  flow  In  the  seal  gap  satisfies  the  same  equations  as  the  flow  In  journal 
bearings.  If  only  one  seal  surface  is  rough »  the  flow  satisfies  (6) 
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where  P  is  the  local  pressure*  x,  y  are  cartesian  coordinates  in  the  plane  of 
the  surface*  h  is  the  variable  spacing  between  the  surfaces*  is  the  viscosity* 
and  U  is  the  local  velocity  of  the  smooth  surface  in  a  coordinate  system  where 
the  rou^  surface  is  at  rest.  The  local  average  fluid  velocities  are 
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We  ask  th«  question. whether  a  steady  vapor  bubble  is  possible  in  the  flow 
described  by  these  equations?  These  equations  do  not  include  details  of  the 
bubble  shape  in  the  third  dimension.  The  effect  of  this  can  be  accounted  for 
to  some  degree  of  approximation  by  imposing  a  pressure  jump  across  the  vapor- 
liquid  interface  equal  to  the  product  of  the  surface  tension  and  the  curvature. 
For  a  steady  vapor  bubble*  the  requirement  at  the  interface  is  that  the  pressure 
is  constant*  and  that  the  interface  is  a  stream  line.  The  bubble  must  also  be 
a  local  minimum  of  pressure*  so  that  the  pressure  must  rise  as  one  moves  away 
from  the  bubble.  This  pressure  gradient  always  drives  the  liquid  toward  the 
vapor  bubble.  The  moving  boundary  always  tries  to  carry  the  liquid  in  the 
direction  of  its  motion.  These  two  forces  can  be  in  balance  for  a  streamline- 


bubble  boundary  only  if  the  normal  direction  of  the  bubble  surface  in  the  direction 


from  the  vapor  to  the  liquid  has  a  component  in  the  direction  of  the  motion 
of  the  moving  boundary*  Ihis  implies  that  the  only  steady  bubbles  which  are 
possible  are  those  whose  upstream  surface  is  a  solid  asperity  that  excludes  the 
fluid  (h*0).  Separation  of  the  flow  into  vapor  and  liquid  layers  in  the  third 
dimension  is  not  likely  in  the  steady  state  unless  the  vapor  in  composed  of  a 
permanent  gas. 

If  the  solution  of  Eq.  (1)  gives  a  pressure  less  than  the  vapor  pressure 
of  the  liquid,  and  the  conditions  stated  above  are  not  met,  the  flow  must 
violate  the  equation  either  by  becoming  unsteady  or  by  developing  a  strong  three 
dimensional  component.  Perhaps  convected  bubbles  appear  which  can  be  carried  to 
high  pressure  regions  where  they  collapse. 

For  flows  where  steady  vapor  bubbles  are  assumed  to  exist,  the  pressure  dis 
tribution  on  the  surface  is  different  than  that  which  would  exist  if  no  bubbles 
were  present.  Calculations  of  the  effect  of  this  steady  cavitation  on  the 
leakage  rate  and  the  forces  on  a  seal  surface  are  presented  in  Section  VII. 


IV.  BUBBLE  DYNAMICS  IN  CONFINED  SPACES 


Some  idea  of  the  possible  pressures  that  can  be  produced  by  bubble  collapse 
In  the  seal  gap  can  be  obtained  by  examining  limiting  forms  of  the  equations  of 
motion  for  bubble  growth.  The  Invlscld  limit  Is  appropriate  when  changes  occur 
so  rapidly  that  viscous  shears  do  not  have  time  to  diffuse  across  the  gap.  In 
this  limit,  the  dynamic  equation  for  a  bubble  In  a  constant  thickness  gap  can  be 
obtained  by  a  simple  extension  of  the  equation  for  a  spherical  bubble  In  an 
Infinite  space  (7) .  For  the  limit  that  the  bubble  radius  Is  much  greater  than 
the  gap  thickness,  the  equation  Is 


Here  R  Is  the  radius  of  the  disk  shaped  bubble,  U  Is  dR/dt,  r*  Is  some  outer 
boxmdary  at  a  distance  much  greater  than  R,  h  Is  the  gap  thickness,  a  Is  the 
surface  tension  constant,  P  Is  the  vapor  pressure  of  the  liquid,  p  Is  the  fluid 

yf  *  CO 

pressure  at  large  distances  from  the  bubble,  p  Is  the  fluid  density,  and  NT  Is 
related  to  the  amount  of  Inert  gas  within  the  bubble.  Values  of  R  and  (P  P  ) 
which  make  the  right  hand  side  equal  to  zero  are  possible  equilibrium  sizes  for 
a  given  nucleus.  Figure  2  shows  this  relation  graphically.  For  a  given  NT, 
there  Is  some  pressure  below  which  no  equilibrium  size  can  exist.  For  values  of 
R,  (P  -  P  )  off  the  equilibrium  curve,  Eq.  (4)  shows  hw  R  changes  with  time. 

xo  V 

If  we  define 
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the  equation  can  be  written 
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(5) 


T  and  R  are  characteristic  length  and  time  scales.  The  solution  of  this 
o  o 

equation  for  two  different  pressure  histories  are  shown  In  Fig.  3  and  Fig.  4. 

In  Fig.  3  the  change  In  R  for  a  step  Increase  In  pressure  Is  shown.  The  Initial 

conditions  were  that  R  Is  an  equilibrium  value  at  t'»  0,  and  dR/dt  >0.  A  step 

Increase  In  pressure  is  applied  and  R  Is  seen  to  oscillate  with  a  period  dependent 

on  the  size  of  the  pressure  jump.  In  Fig.  4,  a  low  pressure  Is  first  applied  so 

that  the  bubble  grows,  and  then  a  high  pressure  Is  applied.  Here,  the  bubble 

collapses  to  zero  size.  This  behavior  Is  similar  to  that  of  a  spherical  bubble 

and  Is  known  to  lead  to  high  pressures  and  cavitation  damage  to  surfaces.  This 

2 

type  of  behavior  arises  from  Inertia  effects  that  are  represented  by  the  U  terms 
In  the  dynamical  equation. 

The  other  limiting  case  for  bubble  growth  corresponds  to  such  slow  changes, 
that  viscous  forces  balance  the  pressure  forces  and  the  Inertia  Is  of  no  impor¬ 
tance. 

For  this  limit,  the  flow  In  the  neighborhood  of  tl. ^  liquid-vapor  interface 
presents  some  difficulties.  A  complex  interaction  of  surface  tension,  viscous 
stresses,  and  pressure  determine  the  Interface  shape,  the  thickness  of  the  fluid 
layer  adhering  to  the  solid  surfaces,  and  the  pressure  changes  In  the  neighborhood 
of  the  bubble  (  8,9).  A  qualitative  view  of  the  flow  in  this  limit  can  be 
obtained  by  balancing  the  bubble  pressure  with  a  quasl-statlc  viscous  flow.  An 
approximate  equation  Is 


This  equation  does  not  have  the  non-linear  terms  that  gave  the  solutions  of  the 
Invlscld  equation  their  violent  behavior  on  collapse.  If  the  flow  In  seal  gaps 
Is  controlled  to  a  good  approximation  by  this  equation,  then  the  rapid  collapse 


of  bubbles  leading  to  cavitation  damage  Is  not  likely.  The  time  scale  for  growth 
and  collapse  becomes  the  crucial  element  In  determining  which  type  of  flow  field 
will  exist,  and  the  physical  basis  for  estimating  this  becomes  the  goal  of  this 
work. 


V.  STRESSES  IN  THE  SOLID 


Cavitation  damage  results  from  the  stresses  applied  to  the  solid  surface  by 
the  flow  field  associated  with  a  collapsing  bubble.  For  extended  flows  only  a 
small  fraction  of  the  collapsing  bubbles  produce  damage  (10).  This  Is  thought 
to  be  due  to  the  requirements  that  damaging  collapse  must  occur  near  the  surface, 
and  also  must  be  of  the  type  (asymmetric)  and  orientation  that  can  cause  damage. 
These  two  requirements  strongly  restrict  the  number  of  damaging  blows. 

Cavitation  In  a  gap  always  occurs  In  close  proximity  to  the  surface  and 
only  one  type  of  collapse  Is  likely.  If  damaging  collapse  occurs  at  all.  It  Is 
likely  that  almost  every  bubble  will  produce  damage. 

The  stresses  In  a  solid  produced  by  an  axlsymmetrlc  distribution  of  pressure 
on  the  surface  can  be  computed  within  the  elastic  limit  by  existing  theory  (11) . 
Knowledge  of  the  pressure  field  generated  by  the  fluid  motion  will  easily  lead 
to  an  assessment  of  Its  damage  producing  potential.  The  shear  stress  on  the 
axis  produced  by  two  different  dlstlrbutions  of  surface  pressure  Is  shown  In 
Fig.  5.  One  Is  for  a  uniform  pressure  within  a  circle,  and  the  other  Is  for  a 
ring  shaped  distribution.  The  latter  distribution  is  more  typical  of  that 
associated  with  bubble  collapse.  In  both  examples,  the  shear  stress  on  the  axis 
Is  maximum  at  a  short  distance  below  the  surface,  and  Is  a  fraction  of  the  sur¬ 
face  pressure.  Other  distributions  of  pressure  give  similar  results. 


VI.  EXFERIMEMTAL  PROGRAM 


An  apparatus  was  constructed  that  could  produce  cavitation  In  a  seal  like 
gap,  and  where  the  bubble  growth  and  collapse  process  could  be  produced  In  a 
controlled  and  observable  manner. 

The  Rig 

A  drawing  of  the  rig  test  area  Is  shown  In  Fig.  6.  A  flat  or  shaped  test  sur¬ 
face  ((4)  In  Fig.  6)  25.4  mm  In  diameter  Is  held  against  a  flat,  transparent  quartz 
disk  ((5)  In  Fig.  6)  by  a  hydraulic  cylinder.  The  test  surface  Is  mounted  so  that 
It  can  remain  parallel  to  the  quartz  surface.  A  load  cell  ((7))  supports  the 
quartz  and  measures  the  force  applied  to  the  surface.  The  force  between  the 
disk  and  the  surface  can  be  varied  In  time  with  a  frequency  up  to  200  Hz  and 
with  an  amplitude  of  the  steady  component  plus  the  oscillating  component  of  up 
to  2000  lbs  (8910N) .  The  test  surface  has  a  central  hole  through  which  liquid 
can  be  forced  with  a  pressure  up  to  2000  psl  (13.8  x  10^  Pa).  A  photograph  of 
the  test  area  Is  shown  In  Fig.  7. 

The  hydraulic  cylinder,  and  load  cell  are  part  of  the  closed  loop  drive 
system  shown  schematically  In  Fig.  8.  A  voltage  Into  the  Input  module  causes  a 
force  to  be  applied  to  the  test  surface  with  the  same  time  variation  as  the 
voltage  and  with  a  level  proportional  to  the  voltage.  An  MTS  "SERVAC"  controller 
Is  used.  An  overall  view  of  the  test  facility  Is  shown  In  Fig.  9. 

Test  Procedure 

The  test  surface  Is  held  so  that  It  remains  parallel  to  the  transparent 
surface.  An  Inverted  microscope  permits  observation  and  photography  of  the  Inter¬ 
face  between  the  disk  and  test  surface  by  looking  through  the  flat  quartz  disk. 
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A  proximity  probe  (Bentley-Nevada) Imbedded  in  the  quartz  senses  the  motion  of 
the  test  surface.  A  pressure  transducer  whose  face  can  be  lapped  flush  with 
the  test  surface  is  being  installed  in  the  test  piece. 

For  normal  operation,  a  test  surface  is  made  with  a  stepped  surface  as 
shown  in  Fig.  10.  The  surface  can  be  examined  with  a  tallysurf  profilometer  to 
verify  its  shape.  By  changing  the  liquid  pressure  and  applied  steady  force, 
the  average  gap  size  can  be  set  over  a  limited  range.  The  pressure  distributions 
associated  with  three  different  gap  sizes  are  also  shown  in  Fig.  10,  each  with 
a  different  shape  corresponding  to  a  different  total  force  related  to  the  pres¬ 
sure  integral.  If  this  steady  force  is  modified  by  an  oscillating  force,  the 
pressure  distribution  and  gap  will  vary  accordingly,  and  the  desired  negative 
pressures  can  be  produced  in  a  region  of  the  gap  for  part  of  the  cycle.  Figure 
11  shows  three  Instantaneous  pressure  distributions  for  different  parts  of  a 
cycle  -  when  the  surfaces  are  approaching  one  another,  when  they  are  separating, 
and  when  they  are  at  rest.  The  velocities  shown  correspond  to  a  30  Hz  oscillation 
with  an  amplitude  of  about  half  the  initial  gap. 

The  negative  pressure  region  occurlng  when  the  plates  are  separating  (or 
when  upper  surface  is  moving  upward)  is  the  region  where  cavitation  is  expected 
to  occur.  The  process  of  bubble  growth  and  collapse  in  this  region  can  be  observed 
and  serve  as  a  basis  for  evaluating  the  various  theoretical  models. 

Results 

The  apparatus  has  gone  through  several  generations  of  change,  and  is  now 
rigid  enough  and  with  sufficient  degrees  of  freedom,  to  enable  it  to  perform  as 
planned.  The  control  of  the  pressure  in  the  gap  by  using  the  force  on  the  surface 
as  the  controlled  variable  has  presented  some  problems,  particularly  if  the  force 
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changes  sign  during  part  of  the  cycle.  Control  by  means  of  local  pressure  or  by 
relative  displacement  of  the  surface  Is  now  being  Investigated  and  should  lead 
to  better  defined  conditions  In  the  gap. 

Some  photographs  of  bubbles  produced  by  the  apparatus  are  shown  In  Fig.  12. 

The  top  picture  shows  several  frames  of  a  motion  picture  take  through  a  micro¬ 
scope  at  1000  frames /seconds  of  bubbles  forming  and  collapsing  In  glycerin.  The 
large  viscosity  of  glycerin  slows  down  Its  motion  so  that  the  history  of  the 
motion  Is  accessible.  The  lower  picture  slows  a  still  picture  of  bubbles  pro¬ 
duced  In  water  (1/1000  second  exposure).  The  small  viscosity  and  large  surface 
tension  of  water  makes  It  more  difficult  to  produce  bubbles  and  also  more  difficult 
to  follow  the  details  of  their  history.  At  this  time,  all  we  can  say  Is  cavita¬ 
tion  bubbles  have  been  produced,  but  Insufficient  observations  have  been  obtained 
with  which  develop  understanding  of  the  cavitation  process. 


VII.  THE  FLOW  FACTOR  AND  LOAD  CAPACITY  FOR  CAVITATION  MODEL 
Introduction 

When  sliding  Is  Introduced  In  a  rough  bearing,  there  Is  a  possibility  of 
cavitation  at  the  trailing  edge  of  the  asperities.  Cavitation  Is  the  result 
of  pressure  reductions  In  the  liquid,  more  precisely.  It  occurs  If  the  pres¬ 
sure  Is  reduced  and  maintained  for  sufficient  duration  below  a  certain 
critical  pressure  which  Is  determined  by  the  physical  properties  and  condi¬ 
tions  of  the  liquid.  The  type  of  cavitation  studied  In  this  section  Is  termed 
steady  cavitation;  It  refers  to  the  time  Independent  situation  In  which  the 
liquid  flow  detaches  from  the  rigid  boundary  of  an  Immersed  body. 

When  a  sliding  Is  Introduced  In  a  rough  bearing,  the  pressure  distribu¬ 
tion  about  each  asperity  Is  expected  to  be  antisymmetric.  However,  due  to  the 
fact  that  the  lubricant  cannot  withstand  large  negative  pressures.  It  can 
rupture  at  the  negative  pressure  zones.  There  are  conflicting  views  on 
whether  or  not  this  cavitation  will  always  occur.  For  examples,  for  some 
geometrical  situations.  Burton  (12)  believes  that  cavitation  voids  could  not 
be  expected  to  nucleate  In  the  short  time  available.  However,  Hamilton  et. 
al.  (13)  relates  the  load  support  mechanism  In  face-type  or  axial  seals  to  the 
cavitation.  They  show  that  the  classical  lubrlctlon  theory  does  not  predict 
the  existence  of  load  support  In  the  case  of  flat,  parallel  surfaces  which  are 
separated  by  a  uniform,  steady  film  of  Newtonian  fluids.  But,  since  high 
positive  pressures  overcome  the  negative  pressures  which  are  truncated  due  to 
the  cavities,  these  types  of  seals  do  produce  a  net  load  carrying  capacity. 

In  this  section.  It  Is  assumed  that  cavitation  occurs  whenever  the  fluid 
pressure  drops  below  the  ambient  pressure.  However,  as  Indicated  above, 
althougli  the  pressure  drops  below  the  critical  pressure,  there  Is  some 
question  as  to  whether  cavitation  always  occurs.  Therefore,  one  should  be 
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aware  of  the  fact  that  this  model  nay  under-estimate  the  shear  flow  factor  in 
case  of  cavitation,  which  tends  to  reduce  shear  flow  in  a  bearing. 

A  typical  contact  area  in  a  rough  bearing  and  the  negative  pressure  zones 
is  shown  in  Figure  13.  As  one  may  observe,  most  of  the  cavitation  occurs  at 
the  trailing  edge  of  the  contacts  and  the  rest  occurs  at  the  trailing  edge  of 
the  asperities. 

To  analyze  the  shear  flow  and  load  capacity  in  case  of  steady  cavitation, 
it  is  assumed  that  these  negative  pressure  zones  are,  in  fact,  isobarlc 
regions,  whose  pressure  is  equal  to  the  vapor  pressure  of  the  lubricant. 

Since  the  vapor  pressure  of  lubricant  is  negligible  when  compared  with  the 
pressure  build-up  in  the  bearing,  it  is  taken  to  be  zero. 

The  shear  flow  factor  is  calculated  substituting  this  modified  pres¬ 
sure  into  Equation  A-20  (Appear  in  Appendix  A).  Finally,  the  load  capacity 
for  the  bearing  is  evaluated  using  truncated  pressure  distribution  as 
indicated  in  the  following  sections. 

Shear  Flow  Factor 

The  shear  flow  factors  that  are  presented  in  this  section  are  evaluated 
using  model  problem  3  as  described  in  Appendix  A.  The  presence  of  the 
cavities  affects  the  shear  flow  in  a  bearing  due  to  the  reduction  of  the  areas 
where  lubricant  flows.  The  existence  of  a  cavity  strongly  depends  on  the 

ambient  pressure  and  sliding  speed  of  the  bearing;  therefore,  the  shear 

* 

flow  becomes  a  function  of  ambient  pressure  and  sliding  speed  as  well 

,  *  — 
as  h/o  and  Y.  However,  since  ^  is  non-dlmenslonallzed  by  U  ,  the  shear  flow 

3  8 

factor  becomes 

♦!  “  ♦!  (Wo,Y,P.) 
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The  effects  produced  by  these  paraneters  are  now  discussed.  In  Figure 

14,  a  typical  pressure  distribution  about  an  asperity  Is  shown.  It  Is  evident 

also  from  the  figure  that  as  the  ambient  pressure  Increases  the  change  of 

negative  pressures  near  asperity  decreases.  Therefore,  one  expects  to  get 

higher  flows  as  the  ambient  pressure  rises.  The  effect  of  ambient  pressure  on 

shear  flow  for  Isotropic  surfaces  Is  given  In  Figure  15.  It  Is  observed  that 

* 

as  the  ambient  pressure  Increases,  the  ^  values,  evaluated  considering  cavl- 

9 

tatlon  effects,  approach  asymptotically  the  shear  flow  factors  without  cavita¬ 
tion.  This  approach  of  becomes  much  faster  for  h/o  -  3  and  h/o  >  0.5  .  It 
Is  slower  for  the  Intermediate  values. 

* 

The  faster  convergence  of  shear  flow  factor  at  high  h/o  can  be 
attributed  to  the  offset  of  the  cavitation  effects  due  to  the  other 
parameters.  At  high  h/o  values,  the  volume  occupied  by  cavitation  becomes 
less  dominant  to  the  main  flow.  (In  contrast,  for  small  values  of  h/o  the 
cavitation  effects  become  more  dominant  because  of  an  Increasing  number  of 
contacts  In  the  bearing  area. ) 

Figure  16  shows  the  change  of  with  ambient  pressure  for  Isotropic, 
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* 

transverse  and  longitudinal  roughness.  It  Is  found  that  the  change  In  ^  Is 

B 

larger  for  transverse  roughness  type,  that  Is,  cavitation  Is  more  pronounced 
for  these  types  of  surfaces.  This  Is  mainly  due  to  the  fact  that  the  areas 
behind  the  asperities  which  are  cavltated  are  larger  than  those  of  Isotropic 
and  longitudinal,  the  fluid  must  travel  a  longer  distance  to  fill  the  cavities 
at  the  trailing  edges  of  the  asperities. 

In  Figure  17,  the  variation  of  shear  flow  factor  with  h/o  for  Isotropic 
surfaces  Is  shown.  It  Is  found  that  shear  flow  factor  for  the  cavitation 
model  behaves  the  same  way  as  the  shear  flow  factor  of  the  model  described  In 
Appendix  A  for  non-cavltatlng  flow.  It  approaches  sero  as  h/o  Increases 
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above  3  due  to  the  vanishing  effect  of  roughness;  as  h/o  decreases  it 
increases  to  a  certain  maximum  and  then  sharply  drops  towards  zero.  The 
reason  for  this  behaviour  is  due  to  the  increased  flow  transport  in  the 
valleys.  The  decrease  is  due  to  the  decrease  in  flow  area  because  of 
contacts,  as  described  in  Appendix  A.  However,  it  is  found  that  the  maximum 
points  shift  towards  high  h/o  values  as  ambient  pressure  is  Increased.  When 
cavitation  occurs,  the  flow  is  reduced  due  to  blockade  by  cavltated  regions 
and  contacts.  It  is  also  observed  in  Figure  15  that  for  values  of  h/o  less 
than  1  cavitation  effects  are  less  pronounced  compared  to  h/o  >1  . 

Therefore,  as  h/o  drops  below  1,  the  decrease  in  flow  is  compensated  bF  the 
decrease  in  cavitation.  This  is  because  at  low  ambient  pressures  the 
cavitation  is  more  effective,  and  conqpensatlon  becomes  more  dominant,  delaying 

the  shear  factor  from  dropping  towards  zero. 

* 

Figure  18  shows  the  change  of  ^  with  h/o  for  transversely  and 

o 

longitudinally  oriented  roughness.  The  surfaces  with  longitudinal  roughness 
result  in  a  lower  shear  flow  than  the  surfaces  with  transverse  roughness. 
However,  it  is  clear  from  the  figure  that  transversely  oriented  surfaces  are 
more  sensible  to  ambient  pressure,  since  more  cavitation  occurs  for  that  case. 

Load  Capacity 

The  classical  lubrication  theory  does  not  predict  any  load  support  for 
two  parallel,  smooth  surfaces  sliding  against  each  other.  For  rough  surfaces, 
the  pressure  distribution  about  each  asperity  is  ant l-ssrmme trie,  and  the 
negative  pressures  cancel  the  positive  pressures.  This  results  in  no  net  load 
support,  even  for  rough  surfaces.  However,  since  is  practice  the  negative 
pressure  regions  may  be  cavltated,  rou^  surfaces  can  create  a  net  load 
capacity. 


17 


The  load  capacity  In  a  bearing  la  given  by 


w 


(8) 


where  p’*  Is  the  truncated  pressure  distribution.  Defining  non-dime ns lonal 
pressure  and  load  as 
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Equation  (8)  becomes 


* 

P  - 
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W  -  P*  did?  (10) 

The  results  calctilated  using  Eq.  (10)  are  shown  In  Figure  19. 

Conclusions 

* 

1.  The  shear  flow  factors  ^  calculted  iislng  a  Model  considering  cavlta- 

8 

tlon  (Model  3  In  Appendix  A)  are  found  to  be  dependent  on  ambient  pressure 

,  * 

as  well  as  surface  pattern  parameter  and  h/o  .  The  ^  Increases 

8 

as  P  Increases,  and  otherwise,  exhibits  the  same  behaviour  as  the  non- 
cavltated  model  with  respect  to  changes  In  h/o  and  y. 

2.  The  model  predicts  a  lower  shear  flow  factor  when  cavitation  Is 
considered.  They  approach  asymptotically  to  the  non-cavltated  values  as 
the  ambient  pressure  rises. 


VIII.  COMPARISON  OF  FLOW  FACTORS  FOR  GAUSSIAN  SURFACES 


BETWEEN  THE  OPEN  AND  CLOSED  BOUNDARY  MODELS 


Introduction 

In  this  section,  the  flow  factors  obtained  through  simulation  are 
presented.  Althou^  one  may  use  any  kind  of  surface  roughness,  It  Is 
l]iq>osslble  to  consider  every  possible  roughness  configuration  within  the  scope 
of  this  study.  Therefore,  only  the  Gaussian  surfaces  as  suggested  by  Patlr 
and  Cheng  (14)  have  been  used  to  facilitate  a  direct  comparison  with  their 
results.  Most  of  the  roughness  parameters  of  a  surface  can  be  obtained  from 
two  statistical  functions:  the  frequency  density  of  roughness  heights  and  the 
auto-correlatlon  function  of  the  surface.  Therefore,  a  convenient  way  to 
characterize  surfaces  would  be  to  choose  specific  functions  to  approximate 
real  surfaces.  The  frequency  density  of  the  surfaces  used  Is  chosen  such  that 
It  can-  be  approximated  by  Gaussian  frequency  density  function,  and  the  auto- 
correlatlon  function  of  the  surfaces  Is  chosen  such  that  It  results  In  linear 
auto-correlatlon  functions  for  the  x  and  y  profiles  which  are  reasonably 
approximation  of  the  engineering  surfaces  (15). 

Since  the  engineering  surfaces  exhibit  directional  properties  due  to 
manufacturing  processes  and/or  runnli^-ln,  these  properties  of  the  surfaces 
are  projected  tqr  Y  >  surface  pattern  parameter,  which  Is  the  ratio  of  x  and  y 
correlation  lengths. . 

It  Is  observed  that  shear  and  pressure  flow  factors  are  strongly  depen¬ 
dent  on  the  exact  topography  of  the  surface.  The  values  may  differ 
considerably  even  though  they  are  evaluated  for  statistically  Identical  sur¬ 
faces.  The  scatter  In  these  values  tends  to  Increase  as  h/o  becomes  less  than 
3.  This  Is  mainly  due  to  Increased  dependence  of  flow  factors  on  h/a  and  the 
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directional  properties  of  the  roughness.  Some  typical  examples  of  scattering 
are  given  In  Table  1.  Since  It  Is  found  that  the  flow  factors  do  vary,  we 
decided  to  evaluate  the  flow  factors  for  several  different  but  statistically 
Identical  surfaces  and  average  them.  The  values  presented  In  this  thesis  are 
the  average  of  20  runs,  which  are  believed  to  give  meaningful  averages. 

When  h/o  becomes  less  than  3  the  asperities  start  Interacting.  This 
regime  Is  called  partial  lubrication  regime.  Many  bearings  encountered  In 
practice  operate  In  this  regime.  Given  the  Importance  of  partial 
lubrications.  It  Is  decided  to  consider  only  the  cases  where  h/o  <3. 

However,  when  h/o  Is  less  than  .5  the  asperity  contact  may  become  so  severe 
that  It  forms  barriers,  and  allows  no  flow  through  the  bearing.  This 
behaviour  of  the  contact  may  result  In  unrealistic  values  for  the  flow 
factors.  Therefore,  for  the  above  range  of  the  values  of  h/o,  the  use  of  this 
model  Is  not  advised. 

Pressure  Flow  Factor 

The  pressure  flow  factors  presented  In  this  sectln  are  obtained  using  the 
model  problem  1  described  In  Appendix  A.  They  are  usually  an  average  of 
twenty  runs  which  use  different  but  statistically  Identical  surfaces  as 
Indicated  In  the  previous  section.  It  Is  evident  that  Is  not  a  function  of 
velocity  since  It  does  not  appear  In  the  formulation  of  model  problem  1,  and 
It  Is  also  not  a  function  of  the  nominal  pressure  gradient  due  to  the  fact 
that  Is  normalized  by  this  quantity  (Eq.  A-16).  The  pressure  flow  factor 
Is  then  only  a  function  of  the  film  thickness  —  standard  deviation  of 
roughness  ratio  and  surface  pattern  parameter  of  the  combined  roughnes,  that 

i*» 
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In  Patlr  and  Cheng's  work  (14),  a  no  flow  boundary  condition  is  used  for 
the  sides  of  the  element.  This  boundary  condition  is  Justified  if  the 
elements  used  in  the  average  Reynolds  Equation  contain  a  large  number  of 
asperities.  However,  it  is  of  Interest  to  determine  the  effect  of  an  open 
type  side  boundary,  as  shown  in  Figure  21,  on  the  pressure  as  well  as  shear 
flow  factors. 

The  results  in  this  section  are  designed  primarily  to  examine  the  differ¬ 
ence  between  the  closed  and  open  boundary  conditions.  They  are  presented  as 
functions  of  h/o  and  y  .  In  addition,  the  effect  of  the  ratio  of  the  element 
length  to  asperity  length  is  studied. 

In  Figure  20,  values  are  plotted  against  the  h/a  for  y  “  1  .  The 
solid  line  in  the  figure  represents  the  average  values,  while  the  dotted  lines 
show  the  minimum  and  maximum  deviation  from  the  mean.  As  it  is  also  evident 
from  the  figure  that  the  scatter  in  these  values  Increases  as  h/a  decreases, 
this  is  mainly  due  to  the  fact  that  a  small  deviation  from  isotropy  results  in 
a  large  change  in  .  Therefore,  the  flow  factors  obtained  from  different 
surfaces  scatter. 

Figure  21  shows  the  average  pressure  flow  factor  values  with  open  side 
boundary  for  isotropic,  transverse  and  longitudinal  roughness.  For  surfaces 
having  larger  correlation  lengths  in  the  longitudinal  direction  (y  >  !)>  the 
pressure  flow  is  enhanced.  Therefore,  is  greater  than  unity,  but  for 
isotropic  and  transverse  roughness,  is  smaller  than  one  and  shows  a 
decreasing  trend  when  h/a  decreases  from  3  to  about  1.  For  h/a  <  1,,  an 
uptrend  for  is  evidenced. 

Figure  22  shows  a  comparison  of  between  the  open  boundary  solution  and 
the  closed  boundary  solution.  It  is  seen  that  the  open  boundary  solution 
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yields  a  higher  value  of  particularly  for  low  values  of  h/a*  This  discre¬ 
pancy  seems  to  pose  a  question  on  whether  the  close  boundary  solution  Is  a 
good  approximation  for  calculation  of  these  flow  factors.  This  question  was 
resolved  some  additional  runs  on  the  effect  of  on  the  pressure  flow 

factor  .  The  ratio  represents  approximately  the  number  of  asperities 

In  the  simulated  region.  It  shows  very  clearly  that  the  open  boundary 
solutions  are  affected  by  the  value  of  whereas  the  close  boundary 

solution  Is  totally  Insensitive  to  any  change  of  ^ 

large,  the  two  solutions  appear  to  converge.  This  result  suggests  that  the 
closed  boundary  solution  Is  Indeed  a  good  approximation  for  calculation  of 
flow  factors  as  long  as  there  are  10-15  or  more  asperities  within  the 
elemental  area  for  which  (p  or  ^  will  be  used.  Recently  Tondor  (16)  and 
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also  Teale  and  Lebeck  (17)  also  used  an  approach  similar  to  Patlr  and  Cheng's 
to  calculate  pressure  flow  factors  for  Isotropic. 

Tondor  obtained  values  of  the  pressure  flow  factors  surfaces  much  higher 
than  those  calculated  by  Patlr  and  Cheng.  Tondor 's  model  corresponding  to  the 
open  boundary  solution  but  with  only  one  or  two  asperities  In  the  region  of 
simulation.  The  results  shown  In  Figure  23  explains  that  the  over  estimate 
of  from  Tondor' s  model  Is  mainly  due  to  the  small  sample  of  asperities  In 
the  simulated  region.  It  Is  believed  that  for  the  average  flow  model  to  be 
meaningful  there  must  be  sufficient  asperities  In  the  simulated  region  to 
account  for  the  Integrated  effect.  For  this  reason,  It  appears  that  the  close 
boundary  solution  should  be  the  most  meaningful  results  for  the  average  flow 
model  Reynolds  equation. 
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Shear  Flow  Factor 

The  shear  flow  factors  obtained  using  the  open-boundary  model  are 
presented  here.  The  values  are  computed  for  twenty  surfaces  and  then 
averaged.  Similar  to  the  pressure  flow  factor,  the  shear  flow  factor  Is  also 
a  function  of  film  thickness  and  the  roughness  parameters  alone.  However, 
unlike  which  only  depends  on  the  statistics  of  the  combined  roughness,  the 
shear  flow  factor  depends  on  the  statistical  parameters  of  rouoghness  of  both 
surfaces  separately.  If  we  consider  a  bearing  In  which  surface  1  Is  moving 
and  surface  2  stationary,  and  If  the  two  surfaces  are  statistically  Identical, 
there  will  be  no  net  flow  due  to  the  combined  effect  of  sliding  and 
roughness.  If  the  rough  surface  Is  moving,  the  fluid  carried  In  the  valleys 
will  result  In  an  additional  flow  transport,  and  yield  a  positive  ^  .  On  the 
other  hand.  If  the  smooth  surface  Is  moving,  the  fluid,  trapped  In  the  valleys 
of  the  rough  surface,  will  reduce, the  flow  and  produce  a  negative  ^  . 

Although  the  shear  flow  factor  depends  on  the  surface  roughness  parameter  of 
both  surfaces,  since.  In  this  study,  the  values  are  obtained  by  using  a  smooth 
surface  and  a  set  of  statistically  Identical  surfaces,  this  dependence  Is  not 
cosldered  as  a  parameter. 

In  Figure  2A,  Is  plotted  against  the  nominal  film  thickness  for 
8 

different  surface  pattern  parameters.  As  expected,  approaches  zero 

s 

as  h/o  Increases  because  of  the  vanishing  effect  of  roughness;  however,  the 
flow  transport  due  to  the  roughness  reaches  a  maximum  somewhere 
around  h/a  *  1  ,  and  then  drops  rapidly  towards  zero.  This  behaviour  Is 
related  to  the  contacts  at  partial  lubrication  regime.  Since  contacts  do  not 
permit  flow,  Is  decreased  due  to  the  decrease  In  available  flow  areas 
as  h/o  drops  below  1. 


since  the  longitudinal  roughness  allows  flow  between  valleys  the  flow 

transported  by  the  valleys  of  the  roughness  for  the  longitudinal  roughness  Is 

less  than  that  transported  by  the  transverse  roughness.  The  ^  curves  for  the 

8 

transverse  roughness  are  at  a  much  higher  level  than  the  <p  curves  for  either 
isotropic  or  longitudinal  roughness.  The  use  of  an  open  boundary  on  the  sides 
produces  shear  flow  factors  slightly  lower  than  those  obtained  with  a  closed 
boundary  condition  (Fig.  25).  As  seen  In  Figure  26,  the  open  boundary 
solution  is  sensitive  to  the  average  number  of  asperities  with  the  simulated 
region,  >  As  this  number  Increases,  there  Is  hardly  any  difference 

between  the  open  and  closed  boundary  solutions. 


IX  SIMMARY 


Work  initiated  in  understanding  the  cavitation  phenomenon  in  a  narrow  gap 
between  sealing  surfaces  has  led  to  the  following  results: 

a)  An  experimental  rig  consisting  of  a  hydrostatic  circular  step  pad  was 
constructed  to  simulate  the  cavitation  phenomenon  under  a  water 
pressure  up  to  2000  psl  comparable  to  those  in  the  submarine  stern- 
tube  seals.  Cavitation  bubbles  were  observed  for  water  glycerine, 
and  their  mixtures. 

b)  Several  analytical  models  were  studied,  and  the  physical  parameters 
influencing  bubble  growth  and  collapse  were  identified.  An  analysis 
was  also  made  to  determine  the  axlsymmetrlc  stress  in  the  solid 
induced  by  a  sudden  pressure  distribution  during  the  collapsing  of  a 
bubble . 

c)  An  existing  thln-fllm  analysis  for  determining  the  flow  factors  used 
in  the  hydrodynamic  analysis  of  rough  surfaces  was  extended  to 

Include  the  steady-state  cavitation  effects.  Results  showed  that 

* 

1)  The  shear  flow  factors  ^  calculated  using  a  model  considering 

s 

cavitation  (Model  3  in  Appendix  A)  are  found  to  be  dependent  on 

ambient  pressure  P  as  well  as  surface  pattern  parameter  y  and  the 

* 

film  thickness  to  roughness  ratio,  h/o.  The  ^  Increase  as  P 

8  Sl 

increases,  and  otherwise,  exhibits  the  same  behaviour  as  the  non- 
cavltated  model  with  respect  to  changes  in  h/o  and  y> 

2)  The  model  predicts  a  lower  shear  flow  factor  when  cavitation  is 
considered.  They  approach  asymptotically  to  the  non-cavltated  values 
as  the  ambient  pressure  rises. 
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TABLE  1 

SCATTER  IN  VALUES 


h/a 

X  -X 

X  y 

M-N 

Y 

6  giax 

X 

average 

min 

X 

.75 

3 

25 

1 

.99253 

.90403 

.85333 

1 

3 

25 

1 

1.03306 

.85140 

.73472 

1.5 

3 

25 

1 

1.11279 

.83785 

.62215 

2 

3 

25 

1 

1.36528 
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.53607 

3 

3 

25 

1 
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1.07882 

.36928 
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Fig.  3.  Bubble  response  to  a  step  Increase  in  pressure.  Initial 

R  ■  1.1m»  'll*  1.25  h.  Fluid  is  water.  Each  curve  corresponds 
to  a  different  pressure  Increment.  l.Ap"  0.1  MPa; 

Ap>  0.5  MPa;  3.Ap  «1.0  MPa. 
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Fig.  6.  Drawing  of  the  test  region  of  the  cavltt.  m  apparatus.  Dimensions 
are  in  inches.  4  is  the  test  surface.  5  is  the  quartz  surface 
mounted  on  the  load  cell. 
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Static  pressure  distribution  for  a  stepped  test  surface, 
force  on  the  test  surface  depends  on  the  gap. 
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Fig.  12.  Micrographs  o£ bubbles  appearing  In  the  gap  when  the  gap  height  changes 
with  time,  a)  shows  a  time  sequence  at  1000  frames/ sec  of  bubbles 
In  a  water-glycerin  solution;  b)  shows  bubbles  In  pure  water. 
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Figure  16.  Change  of  shear  flow  factor  by  ambient  pressure  for  isotropic, 
transverse  and  longitudinal  surfaces. 
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20.  Scatter  in  pressure  flow  factor  obtained  from  different  but 
statistically  Identical  surfaces. 
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Figure  22.  Comparison  of  pressure  flow  factor  with  Patlr-Cheng  Model. 


Figure  23.  Change  of  pressure  flow  factor  with  number  of  asperities 
per  simulation. 
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Figure  25.  Comparison  of  shear  flow  factors  with  closed  boundary  model 
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Figure  26.  Change  of  0^  with  niunber  of  asperities  in  an  element. 


APPENDIX  B 


AVERAGE  REYNOLDS  EQUATION 


For  Reynold's  roughness,  a  term  given  to  roughness  with  small 
slopes,  the  average  Reynold's  equation 

3  JL/lll.  i£'>  ^l''^2  2l!l  +  ^  A 

3x42p  ax''  3y'‘12u  3y''  2  3x  3t 

Is  applicable.  Since  most  of  the  engineering  surfaces  are  of  that  type 
of  roughness,  we  assumed  that  local  pressure  in  a  rough  bearing  is 
governed  by  Equation  A.l. 

The  mean  pressure  is  usually  the  desired  quantity  In  a  bearing. 
However,  since  the  H^  is  a  random  quantity  in  Equation  A.l,  the  local 
pressure  la  also  a  random  quantity.  Therefore,  an  average  Reynolds 
equation  Is  derived  In  terms  of  average  film  thickness  and  average  pres¬ 
sure  which  governs  the  flow  in  a  rough  bearing. 

Before  deriving  the  average  Reynolds  equation,  one  should 
analyze  the  expected  flows. 

The  local  oil  flows  in  a  rough  bearing  are  given  by 


56 


The  film  geometry  randomly  varies  In  a  rough  bearing.  Therefore, 
the  local  oil  flows  are  also  random  quantities.  To  obtain  an  average 
oil  flow,  consider  an  element  with  area  AxAy  (Figure  A.l).  Note  that 
this  element  should  contain  sufficient  numbers  of  asperities  and  still 
should  be  small  when  compared  with  the  dimension  of  bearing.  In  this 
case,  the  expected  flow  has  small  variances  In  this  control  volume. 
Therefore,  they  can  be  obtained  by  averaging  along  the  length  of  the 
element.  That  Is, 

_  -  y+Ay 

_  .  x+Ax 


where  E  Is  the  expectancy  operator. 

In  order  to  relate  these  mean  flows  to  mean  quantities  like 
mean  height  and  mean  pressure,  Patlr-Cheng  defined  empirical  factors 
such  that  expected  unit  flows  are  given  by 


■  12u  3x  ^  2 


2  T-  ^l'^2 
-  h_  +  — 5 —  aij) 


3 


A.  6 


-  e  JlllE 

”^y  12u  3y 


A.  7 


where  p  Is  the  mean  pressure  and  h^  Is  the  mean  gap. 

In  Equation  A. 6  the  flow  In  x  direction  consists  of  three  parts. 

The  first  term  stands  for  the  average  flow  due  to  the  pressure  gradient 

In  x  direction.  The  symbols  ^  and  (p  are  called  pressure  flow  factors, 

X  y 


and  can  be  thought  as  the  correction  factor  between  a  rough  and  smooth 

bearing  which  has  the  same  nominal  geometry. 

The  second  term  In  Equation  A. 6  represents  the  flow  due  to  the 

entrainment  velocity  (U^  +  This  term  Is  absent  In  Equation  A.  7 

because  there  Is  no  velocity  In  this  direction. 

The  third  term  In  Equation  A. 6  arises  due  to  the  combined  effect 

of  roughness  and  sliding.  It  represents  the  additional  flow  transport  In 

the  valleys  of  the  rough  surface.  The  term  compares  the  flow  In  a 

s 

rough  bearing  with  a  smooth  one  In  case  of  sliding.  The  fluid  carried 
In  the  valleys  of  the  moving  rough  surface  helps  to  transport  the  flow 
in  the  gap  between  the  two  surfaces.  On  the  other  hand.  If  a  smooth 
surface  Is  sliding  against  a  stationary  surface,  the  flow  transport  will 
be  obstructed  by  the  stagnant  fluid  In  the  valleys  of  the  stationary 
surface.  For  this  reason,  the  will  be  negative. 

S 

If  we  write  down  the  mean  flow  balance  on  the  control  volume 
we  obtain 


3x  3y  "  ~  3t 


A.8 


That  Is, the  net  flow  leaving  the  volume  Is  equal  to  net  flow  coming  into 
the  volume  and  the  rate  of  change  of  volume. 

Substituting  Equation  A. 6  and  A. 7  into  A.8,  the  following  Is 


obtained: 
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In  order  to  calculate  the  mean  pressure  and  mean  flow,  one  should  com¬ 
pute  the  flow  factors  first,  as  indicated  In  Chapters  I  and  II.  Flow 
factors  have  the  properties 

♦  .  ♦  1  h/o  -*•  » 

*  Y 

(b  •*  0  h/o  « 

s 


Derivation  of  and  through  simulation  (Model  problem  1) : 

In  the  previous  section.  It  is  assumed  that  the  bearing  consists 
of  small,  rectangular  bearings  with  area  6A^  and  a  constant  nominal  film 
thickness  of  h.  The  partitioned  bearing  approaches  the  geometry  of  the 
real  bearing  as  the  5A^  becomes  smaller.  However,  6A^  should  be  large 
enough  to  contain  sufficient  number  of  asperities. 

Fox  each  bearing  pressure  flow  factors  can  be  calculated  by 
applying  an  arbitrary  pressure  flow  on  the  boundaries  of  the  element  and 
solving  pressure  distribution  and  flow,  then  comparing  It  with  the 
rough  bearing  with  the  same  nominal  geometry.  Hence,  evaluating  for 


different  nominal  film  heights  and  can  be  obtained  as  a  function 
of  h. 


To  obtain  *  ,  Patir  and  Cheng  (14)  assumed  u,*  U,*  U  (pure 

X  J. 

rolling)  and  considered  the  following  model  for  simulation 


3x42w  3x^  3y^l2y  3y^  3x  ^  3t 


A.  10 


where  h.^  ■  h  +  +  ^2  (s*®  Figure  A.l), 


with  boundary  conditions  (see  Figure  A. 2) : 


h 


59  - 


1)  P  •  Pa  X  -  0 

2)  P  -  Pb 


3)  p  -  (— )  Pb  ~  Pa  y  “  0  y  -  l,. 

X  X 


4)  no  flow  at  contact  points. 

If  we  consider  righthand  side  of  Equation  A. 10,  since  h  is 
constant  within  the  element  and  h^  -  h  +  6^^  +  62 


3x  3t  "  "  3x  3t 


A.  11 


The  and  62  sre  time  dependent  due  to  the  motion  of  the  surface; 
therefore  one  can  write  them  as 


y)»  i 


1.2 


A.  12 


and 


36. 


3t 


-  u 


!!i 

3x 


i  -  1.2 


A.  13 


Therefore,  the  righthand  side  of  equation  A. 10  becomes 


36. 


36. 


-  -i  -i 

Then,  Equation  10  becomes: 

h  3  3 

J.r_2_  i£i  +  J_ri  l£i  .  0 

3x42y  3x^  ^  3y^l2y  3y^  ° 

The  6^  and  $2  are  randomly  generated  with  known  statistical 
properties,  as  outlined  in  Patir  (15). 


A.  14 


A.  15 


The  pressure  In  the  simulation  element  is  solved  using  Equation 


A.  13  and  then  is  calculated  using  Equation  A.  6  and  Equation  A.  5. 
Since  "  ^2  *  Equation  A.  6  becomes 


-3 


-  -♦x  1^  ^  ^ 


After  cancellation  it  becomes 


A; 


Jl_  i£ 

12]i  3x 


A.  16 


where 


l£. 

3x 


The  calculation  of  follows  the  same  steps.  The  only  dif¬ 
ference  is  the  flow  in  y-dlrection  sees  a  different  surface  pattern 

than  <P  .  That  is,  (h/o) ,  y)  •  <^  (h/a,  1/y) .  Therefore,  one  should 
X  X  y 

expect  to  get  similar  i  and  i  values  for  isotropic  surfaces. 

X  y 

Since  pressure  flow  factors  depend  precisely  on  the  roughness 
of  each  surface,  in  order  to  obtain  a  meaningful  average  one  should 
solve  the  same  problem  for  several  times  and  average  it. 


Derivation  of  (p  through  simulation  (Model  Problem  2) : 
s 

The  shear  flow  factor,  ,  is  obtained  by  eliminating  <p  and 

S  X 

<p^  from  the  equation  A. 8.  Although  there  is  a  possibility  of  cavita¬ 
tion  in  a  bearing  when  some  sliding  is  introduced,  this  model  problem 
does  not  consider  that.  This  effect  will  be  discussed  in  the  following 
section  in  model  problems. 
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To  eliminate  and  <{i  from  equation  A.  8,  the  rolling  velocity 

X  y 

is  taken  to  be  zero  (l.e.,  »  ^2^  ‘  model  problem  for  <|>g  becomes 

pure  sliding  of  two  nominal  surfaces: 


iE]  +  iE,  .  !lll  M7 

3x^l2y  3x^  3y^l2u  3y^  9t 


-  h  +  +  ^2 


U.  -  -U,  »  1/2 
12  s 


The  boundary  conditions  are: 

1)  P  ■  p.  at  X  -  0,  X  -  L 

A  X 

2)  p  ■  y  ■  0,  y  ■  Ly 

3)  No  flow  at  contacts. 

To  obtain  flow  factors,  an  Infinitesimal  displacement  Is  induced 
on  two  surfaces  and  the  pressure  Is  solved. 


Solving  from  Equation  A. 9, 


^x  "  12u  “  L  L  0  0  12u  3x^ 

X  y 


L  L  ^3 


A.  18 


Since  the  mean  pressure  gradient  is  zero,  and  there  Is  no  rolling 
velocity,  this  expression  is  equal  to  the  additional  flow  transport 
due  to  sliding: 

_  U 


Combining  A.  19  and  A.  18,  <|i  becomes 

s 

^  (-  ^  dxdy 


•  1 


A.  20 


Derivation  of  ^  with  effect  of  cavitation  (Model  Problem  3); 
s 

The  shear  flow  factor  is  evaluated  for  cavitation  using  Equa¬ 
tion  A. 17  of  Model  Problem  2  with  different  boundary  conditions.  It 
is  assumed  that  the  lubricant  will  cavitate  the  regions  where  the 
pressure  is  negative  and  the  pressure  will  be  equal  to  cavity  pressure 
which  is  practically  zero.  The  Model  Problem  3  becomes  the  solution 
of  Equation  A. 17  with  boundary  conditions: 

(1)  p  -  p^  X  «  0  ,  *  " 

(2)  p  -  p^  +  C2  y  *  0  ,  ^  “  S 


(3)  p 


P 

0 


if  p  >  0 
if  p  <  0 


0  <  *  < 

0  <  y  <  Ly 


(4)  No  flow  at  contact  points 


APPENDIX  B 


FINITE  DIFFERENCE  FORMULATION 


The  Reynolds  equation  for  Model  Problem  1  is  non-dimens ionalized 


as: 


where 


2 


_L  ru  ^ 

3y  ‘•“t  3yJ 


0 


P 


p~pb 


‘’t 

«T  “T 


Similarly,  the  Reynolds  equation  for  Model  2  becomes: 


L  2 

y 


3y^®T  3y^ 


at* 


where 

*  V 


B.l 


B.2 


B.3 


B.4 


and  the  Reynolds  equation  for  Model  Problem  3  becomes  similar  to  Equa¬ 
tion  B.3  and  B.4.  However,  since  the  pressure  in  Model  Problem  3  is 
sensitive  to  nominal  film  thickness  and  asperity  size  rather  than  stan¬ 
dard  deviation  of  roughness  and  bearing  length  L  ,  the  above  expression 
for  F  is  non-dlmenslonallzed  as 

h^(p-pJ 

^s  “  6mU  L 

S  X 
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where 


L  *  AxA 

X  X 

and 

p  Ji! 

“  NX 

X 

The  boundary  conditions  are  similarly  normalized.  Since  the 
lefthand  side  of  the  Equations  B.l  and  B.3  are  the  same,  two  problems 
can  be  solved  simultaneously  using  the  same  coefficient  matrix. 

The  finite  difference  equations  for  the  two  model  problems  can 
be  written  In  the  form: 


G. .P.  ,  .+E. ,P  ,+D, -P. .+G.  .+E.  .,,P.  ...  „  _ 

Ij  1-1, j  Ij  l,j-l  Ij  Ij  1-1, j  1+1>J  1,J+1  1,J+1  *  B. 


for  1  =  1,2, . . . ,N 


j  »  1,2, ... ,M 


where 


^Ij  "  “t  1-1/2, j 


^Ij  “  "®T  l,j-l/2 


°ij  "  "^^ij  ®i+i,j  hi  ^l,j+l^ 


To  obtain  that  form  of  equation  both  sides  of  the  equation 


have  been  multiplied  by  -(Ax)  and  the  grids  are  chosen  such  that: 


AxL 


that  Is 


Ax  ■  Ay 

The  coefficient  matrix  of  the  Equation  8.6  may  be  stored  as 


symmetric  banded  matrix  with  Mx  MN  dimensions.  The  finite  difference 


equation  in  matrix  form  becomes 


[c]tp]  -  [f]  B.7 

and  can  be  solved  by  Gaussian  elimination. 

Then  flow  and  shear  flow  factors  can  be  obtained  by  numerical 
Integration  of  Equations  A.  16  and  A. 20  of  Appendix  A. 

To  Include  no  flow  boundary  conditions  at  contact  points,  the 
Ep  at  a  half  grid  point  Is  set  equal  to  a  small  e  whenever  Its  value 
Is  negative.  Although  the  e  does  not  affect  the  pressure  distribution 
through  the  bearing.  It  creates  some  pressure  peaks  under  the  contacts 
to  avoid  the  pressure  gradient  created  only  due  to  the  numerical  reason, 
the  flow  Is  set  equal  to  zero  during  the  integration  process  to  calculate 
flow  factors.  Another  numerical  problem  arises  when  a  pressure  point 
surrounded  by  four  contact  points.  To  avoid  this  behaviour,  these 
equations  are  helpful. 

Finally,  one  should  note  that  G  and  E  in  Equation  B.6  are  the 
third  power  of  the  heights  at  half-grid  points.  These  heights  at  half¬ 
grid  points  cannot  be  obtained  by  averaging  the  heights  at  grid  points 
since  it  is  a  random  quantity.  Therefore,  heights  should  be  generated 
at  half-grid  points  as  well. 
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